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We introduce and discuss a scheme for Cooper-pair pumping. The scheme relies on the coher- 
ent transfer of a superposition of charge states across a superconducting island and is realized by 
adiabatic manipulation of magnetic fluxes. Differently from previous implementations, it does not 
require any modulation of electrostatic potentials. We find a peculiar dependence of the pumped 
charge on the superconducting phase bias across the pump and that an arbitrarily large amount of 
charge can be pumped in a single cycle when the phase bias is ir. We explain these features and 
their relation to the adiabatic theorem. 
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I. INTRODUCTION 

A Cooper-pair pumpP is a superconducting device 
that can be used to transport Cooper pairs by manip- 
ulating some of its parameters in a periodic fashion. 
Cooper-pair pumping has recently attracted considerable 
theoreticaPHS] a nd experimentations! attention. 

Part of this attention stems from the geometric prop- 
erties of the parametric cycle used to perform pump- 
ing. These properties, in turn, leave a distinctive finger- 
print in the pumped charge. The link between geometric 
phases and pumped charge has been established in the 
adiabatic limit, where an explicit rela tion connects the 
pumped charge to the Berry phase^l2U, as first shown in 
Ref. [3] and experimentally demonstrated in Ref. [T7] In 
addition, the breakdown of the adiabatic theorem caused 
by Landau-Zener transitions can be detected as a de- 
crease in the pumped charge^. This offers the opportu- 
nity to develop Landau-Zener-Stiickelberg interferometry 
based on geometric phases^. Finally, it has been pro- 
posed to exploit Cooper-pair pumps for the observation 
of nonabelian geometric phaseS^llil. 

Another reason to study Cooper-pair pumps is that 
they are convenient solid-state implementations of a 
driven quantum two-level system. In the presence of 
a dissipative environment, the pumped charge is deter- 
mined by the quasistationary state reached by the system 
and thus is a sensitive probe of decoherence effects. This 
explains why a Cooper-pair pump was chosen as a "case 
in point" in several theoretical works aimed at studying 
the role of dissipation in driven quantum systems, espe- 
cially i n asso ciation with the breakdown of the adiabatic 
theorerrPMHUl. 

Different types of Coop er-pair pumps have been pro- 
posed and realizec P 5 * 16 * 18 l All these devices comprise the 
same building blocks, namely, superconducting islands 
connected to each other and to superconducting leads 
by Josephson junctions. They are intended to be oper- 
ated in a regime where the charging energy of the islands 
is much larger than the Josephson energies of the cou- 



plings. Thus, at the heart of these implementations is the 
"classical" phenomenon of Coulomb blockade. Pumping 
relies on periodic modulation of electrostatic potentials, 
tuned by gate electrodes which act as pistons in pulling 
Cooper pairs onto and off the islands. As a result, the 
system spends most of its time in a state where Cooper 
pairs are localized, except for the short time required 
for their tunnelling between adjacent islands. The fact 
that a Cooper-pair can be delocalized between, say, an 
island and a lead, is of little relevance for a conventional 
Cooper-pair pump. 

In this Article, we undertake a different approach to 
Cooper-pair pumping, that we call "flux pumping" . It 
is based on the adiabatic, coherent transfer of a super- 
position of charge states across a superconducting island, 
realized by magnetic flux control. In contrast to previous 
implementations, no modulation of electrostatic poten- 
tials is required. The device we consider uses the same 
hardware as the Cooper pair sluice introduced in Ref. 4 
yet pumping is achieved in a completely different manner. 
The dependence of the pumped charge on the phase bias 
across the pump reveals unprecedented features. Among 
them, we find that for a particular choice of the system 
parameters, an arbitrarily large charge can be pumped 
per cycle. However, this is by no means inconsistent as 
at the same time the adiabaticity criterion requires the 
pumping cycle to be correspondingly slow. 

The outline of this paper is as follows. In Sec. [TT} we in- 
troduce the Cooper-pair sluice and set up the theoretical 
framework on which our calculations are based. In the 



core Sec. |III[ we describe flux pumping. In Sec. [TV] we 
characterize the breakdown of the adiabatic limit by per- 
forming numerical simulations with a Master-Equation 
approach. Finally, in Sec. |VJ we summarize our results 
and comment on the feasibility of our proposal. 



II. THE COOPER-PAIR SLUICE 

A schematic drawing of the Cooper-pair sluice is shown 
in Fig. [l](a). It is a fully tunable Cooper-pair transistor, 
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FIG. 1. (Color online) Pumping with magnetic fluxes, (a) 
Schematic circuit of a Cooper pair sluice. A superconduct- 
ing island (green) is coupled to superconducting leads by two 
SQUIDs, acting as tunable Josephson junctions of energy Ji, 
J r . A gate capacitively coupled to the island controls its po- 
larization charge n g . The superconducting phase of the two 
leads is held at a fixed difference tp. (b) Representative time 
modulation of J; , J r leading to flux pumping. The gate posi- 
tion is kept fixed throughout the modulation. 



consisting of a small superconducting island connected 
to leads by two superconducting quantum interference 
devices (SQUIDs). The SQUIDs are controlled indepen- 
dently by adjusting the magnetic fluxes threaded by their 
loops, so that they can serve as Josephson junctions of 
tunable energy J;, J r . A gate electrode capacitively cou- 
pled to the island controls its polarization charge n g . We 
will assume a constant superconducting phase difference 
tp across the device. This can be achieved, e.g., by shunt- 
ing the sluice with a large Josephson junction^. 

We use the sluice in the regime where the charging en- 
ergy Eq is large compared to the tunnel energies Ji and 
J r . We describe the dynamics in the basis of eigenstates 
of charge on the island, and restrict the Hilbert space to 
the states |0) and |1) with no and one excess Cooper pair 
on the island, respectively. In this two-level approxima- 
tion, the Hamiltonian is given in matrix form bjP 



H = 



E, 



£ _ 



-5n g f 
ij- sin j 



J. 



cos ■ 



Ec{\ 



ij- sin ^ 
- Sn a ) 2 



(1) 



where J± — \ ( J/ ± J r ), and Sn g 



the offset be- 



tween the gate charge and the degeneracy point. Follow- 
ing Ref. 1231 we further introduce the following quantities: 



El2 = r V Jf +J'r+ 2 JlJr COS if 
( J r — J I if 

7 = arctan tan — 

' \Jr + Jl 2 
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(2a) 
(2b) 
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As we are operating the pump adiabatically, we will 
study the dynamics in the adiabatic basis consisting of 



instantaneous energy eigenstates. They are given by 

1.9) = j= (V^V |0> + e-^y/l + ^\l)) , (3a) 



|e) = ±= (7l+^|0) -e-^TT^Il)) 



(3b) 



In order to properly account for the pumped charge, it 
is essential to consider the corrections to the instanta- 
neous ground state introduced by the time dependence 
of the Hamiltonian ([T]) . These corrections can be treated 
perturbativeljG^H. p or our purpose, it is sufficient to re- 
tain the first-order correction, which gives for the density 
operator of the sluice: 
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We now turn to the calculation of the pumped charge. 
The current operator for the fc-th SQUID (we set 2e = 1 
henceforth) is formally given by Ik = j^d tpk H, where tp k 
is the phase difference across the fc-th SQUID. Explicitly, 
one finds in the {|0) , |1)} basis: 

Ji ( -r-iv/ 2 
n = ■ 



2i 
2i 







-itp/2 








(5) 
(6) 



The expectation value of the current is given by: Ik = 
Tr(plk) = Id.k + Ip.k , where we have singled out a dy- 
namic contribution I d>k = p gg (g\I k \g) + p ee (e\I k \e), and 
a geometric contribution I p ,k = 2^R.(p ge (e\Ik\g))- The dy- 
namic contribution is readily identified with the super- 
current leaking through the SQUIDs in the presence of 
a phase bias. By contrast, the geometric contribution is 
purely due to the pumping cycle. 

The charge transferred through the fc-th SQUID in a 
period is given by Qtr.fc = Jq h{t)dt. For adiabatic evo- 
lution, charge conservation implies Qtr,i = Qtr.r = Qtr- 
In light of the above discussion, it is possible to define 
a dynamic charge Qd and a geometric (pumped) charge 
Q p so that Qn = Qd + Q P - As shown in Ref. Qd and 
Q p are related to the derivative with respect to tp of the 
dynamic phase 0^ and the geometric (Berry) phase Q b , 
accumulated by the instantaneous ground state along a 
pumping cycle. One has, respectively: 

9Qd (7) 



Qd 



Qp 



dtp 
dtp 



(8) 



Experimentally, Qd and Q p can be distinguished as 
they obey different symmetries. In particular upon re- 
versing the direction of the pumping, Qd is not affected, 
while Q p duly changes its sign. We shall henceforth as- 
sume that such a distinction can be made and only be 
concerned with Q p . 
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III. PUMPING WITH FLUXES 

The principle of flux pumping is illustrated in 
Fig. [I] (b), showing the time evolution of the control pa- 
rameters during a pumping cycle. The charge offset 5n g 
(not shown) is adjusted close to the degeneracy point 
and kept fixed throughout the cycle. At the initial time 
t = 0, the island and the leads are decoupled, with both 
Ji and J r set to their minimum value J m i n - In sector 1 
(0 < t < T/3), the coupling to the left lead is turned on 
by maximizing J;. In sector II (T/3 < t < 2T/3), the 
coupling is swapped from the left to the right lead, in 
such a way that the sum J; + J r is kept constant. Fi- 
nally, in sector III (2T/3 < t < T) J r is turned down to 
Jmi-n, bringing the system back to the initial state. 

This mechanism is reminiscent of the Cooper-pair 
shuttle 2 5 , as both devices rely on the maintenance of a 
quantum superposition of different charge states. There 
are, however, two fundamental differences between the 
Cooper-pair shuttle and the flux pump we propose. First, 
the movable grain acting as the "shuttle" is coupled only 
to one lead at a time, while in our case simultaneous cou- 
pling of the island to both leads is essential to achieve 
a nonvanishing pumped charge. Second, the operation 
of the shuttle relies on accurate control of the contact 
time between shuttle and leads and also of the dynami- 
cal phases accumulated by the wavefunctions along each 
"trip" of the shuttle. By contrast, the charge pumped 
by the flux pump does not depend on the speed at which 
the cycle is performed, as long the adiabatic theorem is 
obeyed. This is a consequence of the geometrical nature 
of the pumped charge. 

In the following, we shall analyze how the pumping 
cycle of Fig. |T|(b) determines the adiabatic dynamics of 
the pump, and hence the pumped charge. We first discuss 
the case (p = 0, as it allows for an intuitive explanation. 
For simplicity in deriving analytical expressions, we also 
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FIG. 2. (Color online) Pumped charge in the adiabatic limit, 
(a) Pumped charge Q p versus offset charge Sn g for different 
values of the superconducting phase bias tp. (b) Q p versus <p 
for different values of Sn„. 



integration of I p ,k [see Sec. |IT] , or obtained by virtue of 
^ [see Appendix [A] . In each case, the result is: 



Q p {p = 0] 



sgn(<5n g ) 1 



(9) 



A. The case tp = 



where we have introduced the ratio 



At t = 0, the island is in a defined charge state (0 if 
5n g < 0, 1 if Sn g > 0), and the energy difference between 
charge states is 2Ec8n g . 



As J ; increases so that J\ ~ 
Ec&ng, the ground state evolves into a superposition of 
charge states. As the drive is adiabatic, charge flows 
from the left lead onto the island in order for the system 
to stay in the instantaneous ground state. The charge 
transferred in this case is simply given by |(l|ff(-|-))| 2 — 
|(I|.g(0))| 2 . In sector II, the swapping of the couplings 
does not change the Hamiltonian ([!]). As a result, the 
ground state of the island does not change and no charge 
is transferred. Finally, in sector III, the same amount of 
charge is released to the right lead, as the system has to 
come back (up to a geometric phase) to the initial state. 

The total pumped charge for the case tp = can be in- 
ferred from this heuristic argument, calculated by direct 



2E c 5n g 



(10) 



In the limit r — > ±oo (corresponding to Ec\Sn g \ <C 
J max ), the absolute value of Q p approaches a maximum 
of half a Cooper pair. This result is approximately valid 
also for a finite J m i n , as long as J mm "C Ec5n g . 

The dependence of Q p on Sn g exhibits a sawtooth 
behavior, as shown by the solid line in Fig. [5Ja). At 
8n g = 0, the pumped charge changes sign discontinu- 
ously. However, when Sn g — > also the minimum en- 
ergy gap AE min = min tS [ ,T] AE(t) = Sn g E c tends to 
zero. This implies that the adiabatic limit, in which the 
present derivation is valid, is only attained for infinitely 
slow evolution. We will return to this point in Sec. |IV| 
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B. The general case 

When <p 7^ 0, the same calculation leading to Q shows 
that the pumped charge in sectors I and III is the same as 
in the case ip = 0. This is only to be expected: as long as 
the island is only coupled to a single lead, the phase dif- 
ference between the leads cannot play any role. The sit- 
uation is different for sector II, where the coupling swap 
now takes place between two leads at different phases. 
As a result, an adjustment of the superconducting order 
parameter on the island is requiredpS This causes an ad- 
ditional geometric current to flow across the sluice, in a 
direction opposite to that of the pumping. 

In Fig. [2|a) we plot Q p versus Sn g for different val- 
ues of ip. For values of ip in the range of and tt/2, Q p 
simply decreases with respect to the case (p = 0. As <p 
is further increased, however, a new trend emerges: Q p 
changes its sign with respect to the unbiased case, except 
in the vicinity of the degeneracy point. The magnitude 
of the counterflowing Q p can well exceed a Cooper pair. 
Finally, at <p = tt the sign of the pumped charge is op- 
posite to that of the unbiased case for all values of Sn g . 
Furthermore, Q p diverges as 1/Sn g for Sn g — > 0. 

The full dependence of Q p on <p is shown in Fig. [2jb) 
for three selected values of Sn g . The reader may guess 
that the integral of each curve in Fig. [5Ja) vanishes. In- 
deed, using ^ we obtain (Q p ) = ^ // Q p (p)d(p = 
0s(27r) — ©b(0). While for the Cooper-pair sluice Ob = 
<p + e sin p (e is a small correction) pi so that (Q p ) = 1, 
in our case Ob(27t) = ©b(0), so that {Q P ) V = 0. This 
implies that flux pumping can only be observed in the 
presence of a well-defined phase bias, for if ip. randomly 
oscillates in time, then no net charge is transferred on 
average. On the other hand, Q p exhibits some degree 
of robustness against small phase fluctuations. In par- 
ticular, for Sn g <C 1, Q p develops a plateau centered at 
ip = 0. This can be seen in the increased flattening of 
the curves with smaller 5n g in Fig. [2lb) (an analytical 
argument is provided in the Appendix). 

We shall now fix our attention on the case ip = n, 
for which we can present analytical results. Upon direct 
integration of the current operator, we obtain for the 
charge pumped in the sector II 



Q P U) [<p = n] = sgn(6n 9 



2\/l + r 2 



(11) 



This must be added to the contribution ([9| from sectors 
I,III, to give the total pumped charge 



Q, 



,[<p = ir} = ^sgn(Sn g ) 



1 + r 2 



(12) 



An alternative derivation of (fl~2] ) involving the Berry 
phase is shown in Appendix |A| On comparing ([9| and 
(12), we see that the pumping direction for ip — tt is al- 
ways opposite to that for ip — 0. From ( 12 ), it is apparent 
that Q p diverges for <p = tt and r — > oo (or Sn g — > 0). 



Notably, one finds that this divergence is not removed 
even when relaxing the constraint J m j n = 0. 

The present results arc supposed to hold in the adia- 
batic limit, which requires \p ge \ -C 1. For ip — tt, \p ge \ is 
maximum at t = T/2, where it attains the value 



3hr 2 

max \p nP \ = 



(13) 



We thus see that \p ge \ also diverges for r — > oo. This 
implies that when approaching the degeneracy point, the 
pumping period should be increased according to T oc r 2 
in order to stay in the adiabatic limit. In other words, 
the increase in Q p comes at the cost of an increasingly 
long T. One can further check that the pumping current 
I p does not diverge at any time. In fact, I p never exceeds 
J max R3 ^jff^pge -C ^V^, that is, much less than the 
critical current of the SQUIDs. 



ADIABATIC BREAKDOWN AND 
DECOHERENCE 



IV. 



In Sec. |III| we carried out our calculations assuming 
adiabatic evolution during the pumping cycle. In this 
limit, the pumped charge is a property of the loop de- 
scribed in the parameter space of the Hamiltonian ([!]). 
As such, it does not depend on the pumping frequency. 
We also pointed out, however, that the validity of the adi- 
abatic theorem requires the condition \p ge \ <C 1 to hold. 
We then warned the reader that in the limit Sn g — > 0, 
due to the vanishing of the instantaneous energy gap at 
t = (and at t = T/2 when p = tt), this condition re- 
quires the pumping cycle to be infinitely slow. As real 
measurements are always performed at finite frequencies, 
this implies that none of the traces of Fig. [5Ja) can be 
exactly reproduced in an experiment in the vicinity of 
Sn g — 0. 

In this Section, we investigate how the breakdown of 
the adiabatic limit affects the pumped charge. As soon 
as we allow nonadiabatic transitions to take place, the 
dynamics of the system b ecomes highly nontrivial. As 
argued elsewher e! 12 * 14 * 22 * ^ , the pumped charge is a sen- 
sitive probe of this dynamics. Indeed, it was exploited 
in Ref. [19] to characterize Landau-Zener transitions in 
the Cooper-pair sluice. Furthermore, the final state of 
the pump at the end of the pumping cycle is in gen- 
eral different from its initial state. As measurements are 
typically averaged over very many cycles, the quantity 
of experimental interest becomes the stationary pumped 
charge Qp*™ determined by the interplay between the 
nonadiabatic drive and the electromagnetic environment 
in which the pump is embedded. 

The numerical results we present are obtained using 
the Master-Equation approach developed in Refs.l2"2"1and 
[23l which consistently accounts for the combined action 
of a quasi-adiabatic drive and decoherence on an open 
quantum two-level system. This formalism is not in- 
tended to address the fully nonadiabatic case; yet, it can 
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FIG. 3. (Color online) Breakdown of the adiabatic limit. Q p 
versus Sn g for <p = (a) and <p = tt (b), for different pumping 
frequencies / (solid lines). The results are obtained by nu- 
merical integration of the Master Equation of Ref. [22] The 
asy mpto tic (adiabatic- limit) Q p , calculated according to Q 
and (12 l, is also shown for comparison (dashed lines). Rele- 
vant parameter values for the pump are: Ec = 1 K, J ma x = 
Q.lEc, Jmin = 0.03J max . For the fictitious environment (see 
Ref.[23]for details): g = 0.02, R = 300 Ml, T = 0, T = 0.4 K. 



be conveniently used to investigate the parameter region 
where the adiabatic limit breaks down. Decoherence (de- 
phasing and relaxation) is modeled by attaching a ficti- 
tious environment to the pump, in the form of a resistor 
capacitively coupled to the central island. This mimics 
the effect of charge noise, which is known to be the first 
cause of decoherence in charge-based devices^. 

A full characterization of the effects of decoherence on 
flux pumping is beyond the scope of this work. Still, 
including decoherence in the model is essential in order 
to reach a quasistationary state. In the following, we 
will consider a zero-temperature enviroment and set the 
coupling parameter to a value small enough that the ef- 
fects of nonadiabatic transitions are not washed out by 
environment-induced relaxation 22 . 

In Fig. pi we plot Qp versus 5n g for the emblematic 
cases tp =~~0 (a) and tp = tt (b). We choose the realis- 
tic device parameters Ec = 1 K, J max = O.lEc, Jmin = 
0.03J max , and explore different pumping frequencies 
(solid lines). The adiabatic-limit predictions for the two 



cases (Eqs. [9] and 12 respectively) are also plotted for 
comparison (dashed lines); notice, however, that they 
were derived in the limit J m ; n — > 0. 

In general, nonadiabatic transitions result in a decrease 
of Q p . This is qualitatively accounted for by the fact that 
the adiabatic excited state of the sluice carries an oppo- 
site Q p with respect to the ground state. For the case 
tp = [Fig. [3][a)], this results in a smearing of the adia- 
batic sawtooth. Note that since we have considered the 
realistic case J m i n 7^ 0, in the limit Sn g — > one still has 
A^min = Jmin, so that the residual coupling partly holds 
back nonadiabatic transitions. Analogous considerations 
can be made for the case <p — n [Fig. [3](b)] . As Sn g is 
reduced, however, the nonadiabatic behavior is no longer 
mitigated by the presence of a finite J m i n (this relates to 
the fact that £ , 12 vanishes for J; = J m when tp = 7T, see 



Eq. 2a). The effect is thus more dramatic, with higher 
frequencies hitting the nonadibatic onset first. We are 
not running the simulation for parameter values too close 
to 5n g — 0, as they fall outside the range of validity of 
our Master Equation. 



V. CONCLUSIONS 

We have presented a scheme for Cooper-pair pumping 
which uses neither a bias voltage nor a modulation of gate 
voltages. This scheme involves the adiabatic and coher- 
ent transfer of a superposition of charge states across a 
superconducting island. It requires a definite phase bias 
to work properly, and presents a peculiar dependence on 
this bias. In particular, when the phase bias is 7r, the 
pumped charge per cycle can become arbitrarily large 
near the gate degeneracy point. The pumped current, 
however, is bounded by the adiabaticity criterion, which 
dictates that in the same regime the pumping frequency 
has to be correspondingly small. 

The realization of this scheme looks feasible, especially 
in li ght of recent results obtained with the Cooper-pair 
sluice^^. To this end, particular care must be taken in 
controlling the phase bias. This can be achieved, e. g., 
by embedding the pump in a fully superconducting loop 
together with a large Josephson junction. In this config- 
uration, the latter would also serve as a current threshold 
detector, as in Ref. [TTJ 

An apparent matter of concern is the fact that at finite 
phase, the supercurrent leaking through the device may 
well exceed the pumped current. For instance, let us take 
the device parameters of Fig. [3] and the pumping cycle of 
Fig. [lja). The mean dynamic current at ip — n/2 can be 
approximated by Id ~ jff,Jmax ~ 350 pA, independent 



of frequency. This is to be compared to the pumping 
current I p = Q p f. At a typical / = 120 MHz, Q p = e 
corresponds to I p f=a 20 pA, so that the pumped current 
accounts for less than 10% of the signal. This is not an 
issue, however, as the supercurrent term is even with re- 
spect to time-reversal symmetry, while the pumped cur- 
rent is odd. As a result, we can determine I p by simply 



G 



subtracting the measured currents when pumping in op- 
posite directions (as done in Ref. fT7|) . 

Still, detecting such a small current circulating in a 
loop may challenge customary techniques. In the search 
for signatures of flux pumping, an important role is likely 
to be played by its distinctive symmetries. Besides time- 
reversal, I p is also an odd function of the gate position 
with respect to degeneracy. Finally, it should not de- 
pend on the direction of the circulating currents in the 
SQUIDs. Altogether, these symmetries may be used to 
rule out the contribution of undcsircd rectification effects, 
possibly originating from spurious inductive or capacitive 
couplings. 
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Appendix A: Flux pumping and Berry phase 

The charge pumped by a superconducting pump in the 
adiabatic limit is linked to the Berry phase b accumu- 
lated by the instantaneous ground state along the pump- 
ing cycle^, as prescribed by For a two-level system 
parametrically driven in closed loop, Qb is proportional 
to the solid angle spanned by the Bloch vector, which 
performs an adiabatic rotation on the Bloch sphere. The 
path drawn by the Bloch vector is shown in Fig. |4^a) for 
the pumping cycle of Fig. [ljb) and a few selected values 
of if. The resulting Qb is plotted versus ip and 5n g in 
Fig. W[b). According to pj), the pumped charge Q p is 
givenby the local slope of the surface plot along the ip 



axis. 



Using (|2j, (|3|, and the definitional, we find 



i j dt{g\g) 
lo 



dt (1 + r)) 7 



(Al) 



For the given pumping cycle, 7 (£) = —p/2 in sector I 
and "f(t) = <p/2 in sector III. The time derivative of 7 
vanishes in these regions. The sudden change of 7 from 
ip/2 to —f/2 at times 0, T, . . . yields a delta function in 
the time derivative, but at that time 1+rj = 0, so that the 



integrand in ( AI I vanishes as well. As a result, the only 



contribution to Qb comes from sector II. The fact that 
sectors I and III do not contribute to Qb is a consequence 
of our choice of adiabatic basis (Eqs. [3j), and is not in 
contrast with the fact that there is a charge flow in sectors 
I, III. Indeed, 0s is only defined for closed loops. It is 
possible to give a gauge-invariant generalization of the 
Berry phase for open loop^3, but we will not need it 
here. 



FIG. 4. (Color online) (a) Plots of the path drawn by the 
ground state of the sluice on the Bloch sphere along a pumping 
cycle, for ip — a (blue), ip = n/2 (red), and ip = tt — a 
(green), with a < 1. The Berry phase is proportional to 
the solid angle spanned by the paths, (b) Berry phase Ob 
versus tp and Sn g . According to (JsJ) , The pumped charge Q p 
is proportional to the slope of the surface plot along the ip 



We will now explicitly calculate Qb and Q p in two 
important cases. 



1. Case <p — 

For definiteness, we assume Sn g < 0. To calculate Q p 
for ip = 0, it is sufficient to expand Qb to first order in 
(p. Up to this order, r\ is time-independent in sector II, 
so that 9 B w |(1 + v)h{2T/3) - 7 (T/3)] = |(1 + 17) V . 
The pumped charge is thus 



Uv = o] = l^-v) 



(A2) 



as we found in 

It is worth noting that near the degeneracy point 



In fact, for sufficiently sma 



<C 1, the validity of ( A2 ) extends to all phases ip =^ n. 



1 Sr 



91 V 



for all times in 



sector II. As a result, the pumped charge is half a Cooper 
pair, as predicted by (A2|. 
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2. Case tp — n 



to give 



Using (|S|) and ( Al I, we find for the pumped charge 

Sn„ 



Q P [<p = tt] = - I 
& J j 



2T/3 

dt 

T/3 



1 



d 2 

x - ; - arctan 



dt dip 



y/Sn* + U - Jr)/{2E C )] 



Using J r + Ji = J max and J r — Ji = 6J roax u/T with 
U = t — T/2, the integral can be evaluated analytically 



T f 

Qp[v = *] = uJ 



T f T ' 6 du 



1 - 



VI + (&ru/T)< 



(A4) 



Vl + r 2 - 1 
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